Intermittent Granular Flow and Clogging with Internal Avalanches 
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The dynamics of intermittent granular flow through an orifice in a granular bin and the associated 
clogging due to formation of arches blocking the outlet, is studied numerically in two-dimensions. 
Our numerical results indicate that for small hole sizes, the system self-organizes to a steady state 
and the distribution of the grain displacements decays as power laws. On the other hand, for large 
holes, the outflow distributions are also observed to follow power law distributions. 

PACS number(s): 05.70.Jk, 64.60 Lx, 74.80.Bj, 46.10. +z 



How much granular mass is expected to flow out 
through a hole at the bottom of a granular bin? If 
the hole is sufficiently small, the outflow of grains stops 
soon due to the formation of a stable arch of grains clog- 
ging the hole. If the arch is broken, the grains in the 
arch become unstable and a cascade of grain displace- 
ments propagates within the system called the 'internal 
avalanche'. As a result, a fresh flow starts through the 
hole, which also eventually stops due to the formation of 
another arch. It is known that the amount of out flowing 
granular mass between successive clogging events varies 
over a wide range, however, on average it depends sen- 
sitively on the hole diameter as compared to the grain 
size. 

In this paper, we study this phenomenon of clogged 
outflow distribution from a granular bin in two- 
dimensions using computer simulations. We observe the 
power law dependence of the distribution. 

This phenomenon is important in the context of recent 
experimental and theoretical studies of the avalanche size 
distributions from the surfaces of a sandpile on an open 
base Research in this direction is motivated by 

the expectation that a sandpile is a simple example of a 
Self-Organized Critical (SOC) system [|[j3). According 
to the ideas of SOC, long range spatio-temporal correla- 
tions may spontaneously emerge in slowly driven dynam- 
ical systems without fine tuning of any control parameter 

Laboratory experiments on sandpiles, however, have 
not always found evidence of criticality in sandpiles. The 
fluctuating outflow of granular mass from a slowly rotat- 
ing horizontal semicircular drum did not show a power 
law decay of its power spectrum as expected from the 
SOC theory M. In a second experiment, the avalanche 
size was directly measured by the outflow amount dur- 
ing the dropping of sand on a sandpile situated on a 
horizontal base. It was observed that the distribution 
of the avalanche size, obeys a scaling behaviour only for 
small piles but not for the large ones ||. However, it 
was claimed that a stretched exponential distribution for 
the avalanches fits the whole range of system sizes ||. 
The reasons for not observing scaling are the existence of 



two angles of repose as observed in Q and also the effect 
of inertia of the sand mass while moving down during 
avalanches. This effect was suppressed in an experiment 
with anisotropic rice grains and scaling was successfully 
observed 0. 




FIG. 1. A granular system of 2000 grains of radii l/\/2 
unit in a bin of width 50 units. The bin has a hole of width 3 
units at the centre. The figure shows the steady state of the 
system after 3000 avalanches. 

In contrast, we believe, internal avalanches are better 
candidates for observing SOC. Due to the high compact- 
ness of grains in the granular material in a bin, a grain 
never gets sufficient time to accelerate much and the ef- 
fect of inertia should be small. Internal avalanches were 
first studied in a block model of granular system by re- 
moving grains situated at the bottom of the bin one after 
another Jq] . While this model and also other cellular au- 
tomata [9| and 'Tetris-like' models Jl(| did exhibit scal- 
ing behaviour for the internal avalanches, the disc models 
|p|, p|jl2| in continuous space, however, did not show suf- 
ficient evidence of SOC. 

We study a two-dimensional granular system with 
grains kept in a vertical rectangular bin. Initially the 
system has no hole at the bottom and we fill the bin 
up to a certain height. A hole of half- width R is then 
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made at the centre of the bottom. Some grains drop out 
through the hole. This flow is stopped when an 'arch' is 
formed, clogging the flow (Fig. 1). In two-dimensions, 
an arch is a sequence of grains, mutually depending on 
one another, such that their weights are balanced by mu- 
tual reaction forces. These arches are formed due to the 
competition among the grains to occupy the same vacant 
region of space. In the steady state, these arches are ob- 
served everywhere. The arch, covering the hole, is broken 
by deleting randomly one of the two grains at the bot- 
tom of the arch. An avalanche of grain movements starts 
within the system. Those grains which reach the bottom 
line with their centres within the hole are considered to 
flow out of the system and are deleted. A deleted grain is 
immediately replaced again randomly on the top surface. 
We measure the number of grains p that flow out of the 
system during an avalanche between two successive clog- 
ging events and study its distributions for various values 
of the hole sizes R. 

The granular system is represented by N hard mono 
disperse discs of radii r. No two grains are allowed to 
overlap but they can touch each other and one can roll 
over the other if possible. A rectangular area of size 
L x x L y on the x — y plane represents our two dimen- 
sional bin. Periodic boundary conditions are used along 
the x direction and gravity acts along the — y direction. 

A 'pseudo-dynamics' Q method is used to study the 
avalanches. In this method we do not solve the classical 
equations of motion for the grain system as should have 
been done in a Molecular Dynamics method. Here, the 
direction of gravity and the local geometrical constraints 
due to the presence of other grains govern the movement 
of a grain. To justify the use of the pseudo-dynamics 
we argue that due to the high compactness of the sys- 
tem a grain never gets sufficient time to accelerate to 
large velocities. During the avalanches, the grains move 
in parallel. This is mimicked in the simulation by dis- 
cretizing the time. In one time unit, all grains are given 
one chance each to move, by selecting them one after an- 
other in any random sequence. A grain can move only 
in two ways: Either it can fall a certain distance or roll 
a certain angle over another grain. We used only one 
parameter 6 to characterize these two cases. Therefore 
6 is the maximum possible free fall a grain can have in 
an unit time. While rolling, the centre of a grain moves 
along the arc of a circle of radius 2r on another grain in 
contact, the maximum length of the arc is also limited to 
S. Therefore, the maximum angle a grain can roll freely 
is: 9 = 5/2r over another grain in contact. We believe 
this dynamics models an assembly of very light grains in 
0+ gravity. 

The presence of other grains in the neighbourhood of a 
particular grain imposes severe restriction to its possible 
movements. This is taken care off in the simulation using 
a search algorithm which digitizes the bin into a square 
grid and associates the serial numbers of grains to the 



primitive cells of the grid. A cell contains at most one 
grain when r > 1/ \/2 is chosen. Then it needs to search 
only 24 nearest neighbouring cells for possible contact 
grains. A grain n in the stable position is supported by 
two other grains with serial numbers and tir on the 
left (L) and right (R). The position of the grain n is then 
updated according to the following prescription: 



• If til = tir = 0, it falls 



• if ri£ ^ but ur — 0, it rolls over til, 

• if til = but riR 0, it rolls over ur, 

• if til 7^ and ur ^ it is stable. 

When the grain n with the centre at (x n , y n ) is allowed 
to fall, it may come in contact with a grain m at (x m , y m ) 
within the distance S. Therefore, during the fall, the co- 
ordinates are updated as: 



z' n = x n and y' n = y m + ^ Ar 2 - (x n - x m ) 2 



if 



Vn ~ y n 



< S otherwise, 
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Similarly the grain n may roll an angle 9' over a grain t 
and come in contact with another grain m at (x m , y m ) . 
The coordinates of the grain n in the new position where 
it touches both t and m are: 
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(x m + x t ) ± (y m - yt)\ 



Ur 2 



Vn 



7j + Vt) T [xn 



x t ) 



Ur 2 



Here d m t is the distance between the centres of the grains 
m and t and the ± signs are for the left and right rolls. 
This position is accepted only if 0' < 6 otherwise the 
grain n rolls down a maximum angle 9 over the grain t. 

The initial configuration of grains within the bin is 
done by releasing grains sequentially one after another 
along vertical trajectories with randomly chosen horizon- 
tal positions. They fall till they touch the growing heap 
and are then allowed to roll down on the surface to their 
local stable positions. This is called the 'ballistic depo- 
sition with re-structuring method' (BDRM) Q. The 
initial pattern does not have arches since a grain while 
rolling down along the surface does not need to compete 
with any other grain. However, the arches are observed 
when the system is disturbed by opening the outlet and 
allowing the grains to flow out. 

After a large number of out-flows, the system reaches 
a steady state. One can characterize this state in many 
ways. We characterize the steady state by the distribu- 
tion of the angles of the contact vectors. When two grains 
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are in contact, we draw two vectors from their centres to 
the contact points. The angle <f> is the angle between the 
contact vector and the +x direction. We measured the 
probability distributions P(<p) in the range 0° — 360° for 
the initial state as well as in the steady state. To arrive 
at the steady state we discard 2N outflows, though the 
transition takes place even earlier than that. 
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FIG. 2. Distribution P(4>) of the angles (f> of the con- 
tact vectors. The plot for the initial distribution is shown 
by black dots. It has four high peaks at 60°, 120°, 240°, 300°. 
The curve with opaque circles corresponds to the steady state. 
The randomness introduced into the system is reflected by the 
drastic reduction of the height of the peaks as well as a broad- 
ening. 

The contact angle distributions are shown in Fig. 2. 
The plot of P{4>) vs. <j) f° r the initial grain distribu- 
tion shows four sharp peaks at 60°, 120°, 240°, 300°. We 
explain them in the following way. If the grains were 
placed at the bottom in a complete orderly way by mutu- 
ally touching one another without any gap, the structure 
generated by the BDRM method would be the hexago- 
nal close packing (HCP) structure where every grain has 
six other grains in contact at the interval of 60°s. How- 
ever, in our case, though we dropped grains randomly at 
the bottom level, the deterministic piling process retains 
some effect of the HCP structure. On the other hand, 
it is also known that the probability of a grain having 
the number of contact neighbours different from four, de- 
creases with the height as a power law |t| , which is cer- 
tainly consistent with the four peaks we obtained. There 
are also two small peaks at 0° and at 180° due to the hor- 
izontal contacts of grains near the base level. This dis- 
tribution, however, changes drastically by reducing the 
heights of the peaks and by broadening their widths in 
the steady state. The randomness introduced into the 



system by randomly selecting grains at the bottom in- 
creases the possibilities of other contact angles also, but 
is not fully capable to totally randomize it to a flat distri- 
bution. We check that after a large number of avalanches, 
this steady state distribution remains unchanged. 




FIG. 3. Distribution P(As) of the displacement As of the 
grains in a granular system in the limit of R — * 0. Curves for 
different sizes with grain numbers N = 300, 625, 2500 and 
10000 are plotted from bottom to top on the left. 



Next we study the statistics of the displacement dis- 
tributions in the granular heap in the limit of R — ► r. In 
this limit, no grain can drop out through the hole. The 
arch is broken by deleting only one grain. The avalanche 
is allowed to propagate within the system. The displace- 
ment As is the absolute magnitude of the displacement of 
a grain before and after an avalanche. We observe that 
there is a huge variation the displacements As of the 
grains. Most of the grains displace very little, whereas 
others have much bigger displacements, but their num- 
bers are small. The displacement distribution P(As) 
is measured in the steady state over a large number of 
avalanches. The lower cut-off of the distribution turned 
out to be strongly dependent on the tolerance factor e 
used for the simulation. If the centres of two grains are 
separated by a distance within r — e and r + e then they 
are considered to be in contact. The upper cut-off of the 
distribution is of the order of unity since when a grain 
is deleted at the bottom, its neighbouring grains drops a 
distance of the order of the grain diameter. 

Systems of four different sizes have been simulated with 
N = 300, 625, 2500 and 10000 in bins with base sizes 
L x = 20, 40, 80 and 160 units respectively. First 57V 
avalanches are thrown away to allow the system to reach 
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the steady state. A list is maintained for the coordinates 
of all grains. This list is compared before and after the 
avalanche to measure As. We observe a nice straight line 
curve over many decades of the P(As) vs. As plot on 
a double logarithmic scale, with little curvature at the 
right edge (Fig. 3). However, the slope of this curve is 
found to increase with the system size. Different extrap- 
olations are tried and the best is obtained with L~ x / 2 
which leads to 0.97 ± 0.05 in the N — > oo limit for the 
slope. We get similar behaviours for the distributions of 
the x and y components of the displacement vectors. 
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FIG. 4. The distribution P(p) of the number of grains p 
that flows out of the system in an internal avalanche between 
two successive clogged situations plotted on a double loga- 
rithmic scale. A bin of size Lx = 50, containing 4000 grains 
of radii r = 1/ v2 was simulated with the hole sizes R = 0.75, 
1, 1.25, 1.5, and 2.0 (curves from left to right). 

Finally we study the outflow distribution and clogging. 
We use a bin of dimensions L x = 50 and L y — 200 units 
containing 4000 grains of radii r = 1/ y/2 with a hole of 
half-width R at the centre of the bottom of the bin. In 
a stable state an arch covers the hole and forbids other 
grains to flow out. The two lower most grains of the arch 
on the two sides of the hole are chosen and one of them 
is selected randomly. By deleting this grain, the arch 
is made unstable, which creates an internal avalanche 
of grain movements within the bin. Grains which come 
down and touch the bottom line of the bin with their cen- 
tres within the hole are removed. While removing such 
a grain, the whole avalanche dynamics is frozen and the 
grains is replaced back on the top surface again by the 
BDRM method. 

Two conical shapes are formed. There are some grains 
which are never disturbed by any avalanche and are lo- 
cated on both sides of the hole. These undisturbed grains 



form a cone. The second cone is formed on the upper 
surface, though the grains were dropped on this surface 
along randomly chosen horizontal coordinates with uni- 
form probability. For a fixed size of the base, the angle 
of the cone is found to increase with the average height 
of the granular column. 

The number of grains p that go out of the system 
in an internal avalanche is counted between successive 
clogging events. Holes of five different sizes, i.e., R — 
0.75,1,1.25,1.5,2 are used. We could generate around 
5 x 10 5 such outflows for each hole size and collected the 
distribution data for P(p). We show the plot of these 
curves in Fig. 4 on a double logarithmic scale. The grain 
which is deleted to break the arch is also counted in p. 
It is found increasingly improbable that the outflow will 
consist of only another grain. Consequently P(l) is found 
to decrease monotonically on increasing R. However, for 
large values of p, the curve falls down. 

All curves of P{p) vs. p plots show straight lines for 
large p values, signifying the power law distributions for 
the avalanche sizes. The slopes of these curves a(R) are 
found to depend strongly on the size of the hole, as well as 
the grain size: P{p) ~ p~ CT (- R < r ). We observe that a(R,r) 
is actually a linear function of (R/r) as <j(R,r) = k{^) 
with k=1.6. The average outflow < p(R/r) > diverges 
exponentially as: exp(a(^r — 1)) with a = 0.7. 

To summarize, we have studied a numerical model of 
intermittent granular flow and clogging with associated 
internal avalanches from a two-dimensional bin using a 
pseudo dynamics method. The flow of grains through 
a hole at the bottom of the bin are repeatedly clogged 
due to the formation of stable arches of grains, which 
we break to start fresh flows. Our numerical simulation 
results indicate that the system self-organizes to a state 
different from the initial state. In the limit of R < r, the 
distribution of granular distributions follow power laws. 
On the other hand when R > r, the granular flow distri- 
bution is also observed to follow power law. 
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